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I. INTRODUCTION 

Quantization is a phenomenon that changes our bright 
classical perspective into a bit uncertain and at first sight 
rather nonintuitive picture. This picture, however, is 
more rigorous than the classical one, possesses many fas- 
cinating features and has produced a lot of successful 
predictions. 

The subtle problem of transition from classical to 
quantal attracts attention from the early days of quan- 
tum mechanics. Over the years various techniques and 
methods for solving this puzzle have been invented. Our 
aim is not to trace back the complete (hi) story of the 
milestone ideas in this field (for the review we refer to 
PQ). What we want to do is to give a concise exposi- 
tion of the method we have developed. However, since 
we have generalized the original Feynman's path integral 
approach, we will recapitulate this approach shortly in 
section ITTJ 

The main goal of our paper is to obtain a func- 
tional integral formula for the quantum propagator which 
would not refer to Lagrangian and/or Hamiltonian func- 
tion. Quantum propagator is the probability amplitude 
A(gi,ti | qo,to) for the transition of the system from 
the initial configuration (qo,to) to the final configuration 
(<7i, ti). We will derive a closed expression for this quan- 
tity starting from the given set of classical dynamical 
equations of motion. 

The proposed method uses functional integration in 
the extended phase space, but instead of integration over 
path histories we introduce integration over stringy sur- 
faces. This crucial element of our approach is explained 
in full detail in sections |III| and |IV| and in appendices 



[A"| and [B] We also make sure that whenever the system 
under consideration becomes Lagrangian, the stringy de- 
scription reduces to the standard one with the Feynman 
path integral. 

In sections [V] and |VI| some simple examples are scruti- 
nized. It is well known that there are classical dynamical 
systems that cannot be described within the traditional 
Lagrangian or Hamiltonian framework. The stringy ap- 
proach enables us to quantize them straightforwardly. In 
section ^ the quantization of a weakly non-Lagrangian 
system is performed. The transition amplitude is com- 
puted for a particle in a conservative field, whose motion 
is damped by a general A-power friction force F = —Kq A . 
The stringy results for A = 1 are compared with the re- 
sults obtained in the approaches by Caldirola-Kanai |2J, 
Bateman [3] and Kostin 0]. 

One can argue against the stringy quantization of dis- 
sipative systems that it is unable to describe decoher- 
ence phenomena. This is true, but the same objection 
can be raised against the generally accepted heuristic 
approaches by the authors cited above, as well as those 
by Dekker [5], Razavy |B], Geicke (7] and others, sim- 
ply because their kinematical and dynamical prerequi- 
sites are different from the prerequisites of the particle- 
plus-environment quantum models. However, a possible 
argument for the stringy quantization is that the particle- 
plus-environment models are not able to handle satisfac- 
torily the case with the friction force —nq A for the gen- 
eral power A. Moreover, they describe a rather different 
phenomenon, namely the quantum Brownian motion for 
which the total force equals — nq + stochastic term. 
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Section VI contains the analysis of a curious two- 
dimensional Douglas system [8] which is strongly non- 
Lagrangian, i.e. not derivable from any sort of La- 
grangian or Hamiltonian. This causes a fundamental 
problem for all conventional quantization methods, but 
can be dealt with in a rather transparent way in the 
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stringy approach. 

In section |VII| conclusion, discussion and outlook are 
collected. The section also includes some comments on 
the relation between the stringy quantization and the 
Caldeira-Leggett model of quantum Brownian motion |9J , 
as well as the influence functional technique by Feynman 
and Vernon |10j . 

Some rather technical material is left to appendices. 
Appendix [X] is devoted to the stringy variational prin- 
ciple, which plays an important role in the motivation 
of our approach. In appendix [S] computational details 
concerning the surface functional integral for quantum 
friction force systems are presented. 

Historically the first attempt at quantization based on 
the dynamical equations of motion belongs to Feynman, 
see A similar problem was considered by Wigner, 

Yang & Feldman, Nelson, Okubo and others, see |12j . 
Among the recent investigations in this field let us men- 
tion the work of Lyakhovich & Sharapov 1 1 3] and Gitman 
& Kupriyanov [2]. They consider the same problem as 
us, but their strategy is different. In our opinion, their 
approach fits much better the context of gauge field dy- 
namics. 

As Ludwig Faddeev noted during the Edward Wit- 
ten's talk at the Mathematical Physics Conference: 
From XX To XXI Century [26], quantization is not a 
science, quantization is an art. Let us believe that the 
quantization method proposed here will fit the Ludwig 's 
dictum and will be meaningful enough to be considered 
artistic. 



This determines the phase of N. 

A miraculous consequence of the definition of the tran- 
sition amplitude (JTJ (not an additional requirement!) is 
that it satisfies the evolutionary chain rule, or Chapman- 
Kolmogorov equation, 

A(gi, h | q , t ) = J dq'Aiquh \ q' , t) A(q', t \ q , t ) . (3) 

The infinitesimal version of this formula is the celebrated 
Schrodinger equation. 

Quantum states of the system are described by the 
square integrable functions with the standard Hilbcrt 
space structure. Physical observables are hermitian op- 
erators acting on such functions. Given the state ^o(q) 
at the initial moment to one is able to predict the state 
^i(q) at any later moment t\ according to the formula 

Ml) - *i(«)= / dq' MQ,hW ,t )* (q') ■ (4) 

In what follows it is assumed that the above concept of 
states, observables and quantum evolution is valid for 
the stringy quantization of non-Lagrangian systems as 
well. The only new element is a modified prescription 
for the evolutionary integral kernel A(qi,t\ | qo,to)- The 
same kinematical prerequisites can be found also in other 
phenomenological approaches [2]- (7]. 

III. ONE STEP BEYOND FEYNMAN 



II. FEYNMAN QUANTUM MECHANICS 

According to Feynman |15| , the probability amplitude 
of the transition of the system from the space-time con- 
figuration (goi^o) to another space-time configuration 

(qi,h) is 

A(g 1 ,t 1 |gb,*o) = ^ J[Dj]exp{^ Jpdq-Hdt}. (1) 

i 

Here the integral is taken over all paths (histories) j(t) = 
(q(t),p(t),t) in the extended phase space [27 , satisfying 
the conditions q(to) — qo and q(t\) = q\. 

The preexponential factor 1/N in the expression ([!]) 
serves just the normalization. To fix it properly we im- 
pose two physical conditions on the transition amplitude. 
First we introduce an integral condition that ensures that 
the total probability is conserved, 



dqi A(q 1 ,t 1 \q ,t ) A(q 1 ,ti \q' ,t ) = S(q -q' ) . (2a) 



This specifies the absolute value of N. Then we add a 
constraint that expresses the obvious fact that no evolu- 
tion takes place if the final time t\ approaches the initial 
time to, 

lim A(q 1 ,t 1 \q ,t ) = S(qx - q ) . (2b) 

ti— >to 



A possible step beyond the theory summarized above 
consists in the elimination of the Hamiltonian function 
H from formula ([TJ. The price to be paid is the replace- 
ment of the path integration by the surface functional 
integration. 

Our aim is to construct the amplitude for the transi- 
tion between (qo,to) and (qi,t±) starting from the classi- 
cal equations of motion (and not from the Hamiltonian 
function which provides them) 



dH 

dp a 



m 



Pa 



dH 
dq a 



= F a . 



(5) 



In the first set of equations, m is the mass of the parti- 
cle. We restricted ourselves to the simplest case of one 
particle, although it is trivial to generalize the theory to 
a system with an arbitrary number of particles. Note 
also that if the particle is unconstrained and we make 
use of Cartesian coordinates, the momenta p a defined in 
^ reduce to p a . 

Suppose that there exists a unique classical trajectory 
in the extended phase space 7 c ;(i) = (qci(t),p c i(t),t), 
connecting the points (ftb^o) an d (<Zi,£i)- Then we can 
assign to any other trajectory ^(t) — (q(t),p(t),t), which 
enters the path integral in 0, two auxiliary curves 

\o(s)=(q ,ic(s),t ), where 7r(0) = p c i (to), tt(1) = p(t ); 
X 1 (s) = (q 1 ,(f)(s),t 1 ), where 0(0) = p c i{h), (j>(l) = p{h) . 
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FIG. 1: Schematic picture of two auxiliary curves Ao(s) and 
Ai (s) which connect the classical history 7^ (t) with the given 
history 7(f) in the extended phase space. A-curves are located 
in the n-dimensional subspaces of the extended phase space 
in which the momenta are varying while the coordinates and 
time are kept fixed. The contour 9E = 7 — Ai — jd + Ao forms 
a boundary for plenty of extended phase space surfaces. One 
of them, denoted as E, is drawn in the figure. 



The curves are parameterized by the parameter s S [0, 1] 
and live in the momentum subsectors of the extended 
phase space with fixed (qo,to) and (qi,ti), see Figure [T| 
Using these definitions one can write [28] 

J p a dq a -Hdt = J p a dq a - Hdt + j) p a dq a - Hdt , (6) 



Tel 



where <9E = 7— Ai — j c i + Ao is a contour in the extended 
phase space consisting of four curves 7(i), Jcl(t), Ao(s), 
Ai(s). The first integral on the right hand side is the 
classical action S c i(qi,ti | qo,to) (its equivalent for the 
non-Lagrangian case will be specified later), while the 
second integral can be rewritten as 



p a dq a - Hdt -- 
as s 



dp a A (dq a ■ 



dH 

dp a 



dt 



dH 



dq a A dt , 



where S is a surface spanning the contour <9E, i. e. a 
map from the parametric space (t, s) € [io,ii] x [0,1] to 
the extended phase space, 

E : (t, s) ^ E(t, s) = (q a {t, s),p a (t, s),t(t, s) = t) , 

satisfying 



E(t,0)=7c,(t) 
E(t,l)=7(t) 



and 



E(t , s) = A (s) 
E(ti,a) = Ai(s) ' 



The surface E can be viewed as a worldsheet of a string, 
therefore we will call the quantization method using such 



surfaces "stringy." Partial derivatives of the Hamiltonian 
function entering the integral over E can be eliminated 
with the help of the equations of motion (JHJ. By doing 
so, one transforms (|6j) into the following form: 



/ 



p a dq a -Hdt=S cl (q 1 ,t 1 \q ,t )+ Q, (7) 



where the two-form fi is defined as 

n = dp a A dq a - ( — dp a - F a dq a 
V m 



A dt. 



(8) 



This two-form is an object in the extended phase space 
and its structure can be read out from the underlying 
equations of motion. Note that for non-potential forces 
the expression (p a /m)dp a — F a dq a does not reduce to 
dH, so that the two-form is not closed. This becomes 
essential in the next subsection. 

It is obvious that for a given pair of histories (7,7c;) 
there exist infinitely many E-surfaces such that 7 — 7 C / C 
<9E. All of them form a set which we will call U(^ j7cI ) 
[25] . Since no E is preferred and J s fl is only boundary 
dependent, it is natural to average the exponent of ^ 
over the whole stringy set li( 7 7ci ). After doing so we 
obtain the identity 



ex P {\j Pa dq« - Hdt) = ^- J [DT] cxp {I J n}, 

where is the cardinality of the stringy set U(^. 7ci ), 
N 7 := #U( 7i7( , i ), and [2)E] is a functional integration 
measure specified in Appendix [B] If no topology-related 
problems arise in the extended phase space, the infinite 
constant N 7 is independent of the history 7. Taking all 
this into account we can rewrite as 



A(<ft.,ti \q ,t ) 



N 



[DE]exp{^ / n}, 



(9) 



where the set U 7ci over which the functional integration 
is carried out contains all strings in the extended phase 
space which are anchored to the given classical trajectory 
Some surfaces from U 7cl are depicted in Figure [2] 



1 ci- 
ln 

into the overall preexponential factor 1/N and the path 
integral over 7's was converted into the surface functional 
integral, as promised earlier, using the identities 



kj^J-ii^ ouiioLm Hum *^~f c i CU.G ucpibicu ill riguic 

Q, the undetermined constant 1/N 7 was absorbec 



Jm |[2)E] 



[DE] 

U «,7 ci ) 



[DE] 



(10) 



IV. ONE MORE STEP BEYOND - 
NON-LAGRANGIAN SYSTEMS 

In case the Hamiltonian function H is given, the for- 
mula ([9]) for the propagator is equivalent to the Feynman 
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feed 



FIG. 2: Four elements (two open, one closed and one shrunk 
to 7 C ;) of U Tc!) anchored to the given classical history jd in 
the extended phase space. The two open surfaces belong to 
two different stringy classes U(~ i7 ,), while the closed and the 
shrunk surface are both elements of the same stringy class 
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(7ci.7 c l)' 



As the worldsheet E varies, the front and rear 



a function (i.e. the forces satisfy the Helmholtz condi- 
tion) then we can define a local auxiliary Hamiltonian 
Haux such that fQ = d(p a dq a - H aux dt). Having H aux 
we can define the auxiliary classical action S aux . The al- 
ternative approach which tries to identify S c ; with S aux 
has, however, several substantial disadvantages. First of 
all, the procedure of finding / (and subsequently H aux 
and S an:r ) is highly ambiguous. Furthermore, S aux lacks 
symmetries that were originally present in the equations 
of motion. Because of these findings we do not follow this 
strategy hereinafter. 

From the physical point of view the dynamical equa- 
tions seem to be more fundamental than their com- 
pact but ambiguous precursors, Hamiltonian and/or La- 
grangian function, see [16]. Since ([£]) requires just the 
knowledge of the equations of motion, it determines the 
transition amplitude in a completely new way. 

Our proposal for the transition amplitude has ap- 
peared here out of thin air. Actually, we were rewriting 
([!} in terms of stringy surfaces and then, when realizing 
that the integrated function can be written without any 
reference to the Hamiltonian, we postulated formula ^ 
to be valid in general. However, this is not the whole 
story. The stringy functional quantization can be moti- 
vated also by the stringy variational principle. Having 
the dynamical equations and initial and final endpoints, 
one can form U 7c! and ft. Then using these objects one 
can introduce the stringy action functional 



boundary curves, denoted in the text Ao and Ar, vary as well. 



§ : IL 



S i y §(S) 



Q. 



(11) 



formula ( 1 1 we have started with. The surface functional 
integral ( 9 ) is however completely free of H and requires 
just the knowledge of the classical equations of motion. 
This observation allows us to postulate |9f as a quan- 
tization tool in situations in which one cannot use the 
standard Hamiltonian approach. We just have to re- 
lax the requirement of the closedness of the two-form 
f2, following from the definition of f2 for Hamiltonian 
systems. This relaxation is what is hidden behind the 
slightly provocative phrase one more step beyond in the 
title of this subsection. 

The problematic part here is the definition of the classi- 
cal action S c i(q\,ti |<70)^o) f° r non-Lagrangian systems. 
Later we will see that in some specific situations this 
quantity can be read out from the structure of the sur- 
face integral 



[DS] ex P {^y il] 



Another possibility which comes to mind is to use the 
integrating factor of the generally non-closed two-form 
O. This means that one will look for a function f(q,p, t) 
on the extended phase space such that d(fft) = 0. If 
the structure of the dynamical equations allows for such 



This is a variational problem with varying boundaries, 
therefore the total variation has two terms. First one 
specifies the boundary and determines the initial equa- 
tions of motion for j c i . The second one specifies the bulk 
of the stationary world-sheet S (which turns out to be 
shrunk into 7^ itself). Moreover, one immediately veri- 
fies that in the special case when = d(p a dq a — Hdt). 
the stringy variational principle reduces (up to an ad- 
ditive constant) to the celebrated Hamilton least action 
principle. Of course, many subtleties were omitted here, 
but all of them can be found in |17J . or in Appendix |A| 
Thus, our variational principle enables us to perform the 
limit h — > 0, in which we recover the original classical 
dynamics as required. 



V. QUANTIZATION OF FRICTION FORCE 
SYSTEMS 

To examine the functionality of the proposed quanti- 
zation method let us first analyze the simplest friction 
force system. It consists of a particle with unit mass, 
moving in one dimension under the action of the conser- 



vative force F 
EI]. Thus, 



-dU/dq and the friction force — Kq J 



-nq' 



'+F n = d(pdq-±p 2 dt-Udt)-Kp A dqAdt. 
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In this example the surface functional integral can be cal- 
culated explicitly (for more detail see |17| . or Appendix 
[B| . In the course of calculation, the surface functional 
integral in the extended phase space reduces to path in- 
tegral in the configuration space, 



[2)S]exp{^y n} oc 

tl 

OC exp {-- J {\qli ~ U{q cl ) - nq c lPci) dt \ x 

to 

ti 

x I [Dg]exp{^ J {\q 2 -U{q)-n qP ^)dt). (12) 

to 



This combined with the formula ([9| suggests that it will 
be convenient to define the classical action as [3T] 

ti 

SAc/(9i>*il*),*o) = y (|<?cz - U{q c i) - Kq c .iPci)dt. (13) 

to 

When doing so, we obtain a controlled cancelation of 
exp{^S^ c i} with the preexponential factor arising from 
(12 1. The final probability amplitude then assumes a 



compact and reasonable form: 
A Astr (qi,ti\q 0l t ) = 

ti 

= — ^ / m exp{ l -j (\q 2 - U(q) - Kqp£)dt} . (14) 

to 

The preexponential factor lfN Astr can in principle be 
obtained by subjecting A.Astr(qi,ti\qoito) to the condi- 
tions 

Let us explain why one should consider the propaga- 
tor formula (14 1 reasonable. In the path integral above 



there appears an unconventional "external source" term 
— Kqptf. Its appearance guarantees that the quantum dy- 
namics governed by A.Astr{qi, t\ \qo, to) transforms into 
classical mechanics in the limit h — > 0. This follows from 
the simple fact that the unique solution of the saddle 



point equation for the path integral ( 14 1 



dU 
dq 



which satisfies the given initial and final conditions 
q(to) — qo an d q(ti) = q%, is the classical trajectory q c i{t). 

The external source term in (14 1 breaks the valid- 
ity of the Chapman-Kolmogorov (memoryless) equation 
|32| . From the physical point of view it is a desired 
phenomenon. The microscopic origin of friction is some 
environmental interaction. This, however, was not ac- 
counted for here explicitly. What has been considered is 
some effective (macroscopic, phenomenological) interac- 
tion emerging on the classical level only. Microscopically 



the system is a part of a larger system and hence it should 
be affected by the memory effect. 

It is not difficult to compute A^ s tr(<Zi> ^i|90j ^o) for 
the general power A and the potentials U(q) — and 
U(q) = \ijj 2 q 2 . In both cases we can carry out the path 
integration in ( 14 ) explicitly to obtain 



A.Astr(qi,ti\qo,t ) = 



— - — exp { ^SAei(qi,h\qo,to) 



(15) 



where the non-Lagrangian action SAci(qi, t\\qo, ^o) is 
given by the expression (13 1. An open problem is to de- 
termine the preexponential factor 1/N^ str (ti — i , g 1; q ) 
for all powers A except for A = 1 (at least the author is 
incapable to do that). If A ^ 1, the preexponential factor 
is apparently dependent not only on the time difference 
ti — to, but also on the endpoint positions qo and q%. 
This hypothesis is supported by the results obtained in 
[T5] , where the path integral for a non-conservative force 
quadratic in velocity is calculated. The normalization 
prefactor which appears there is explicitly dependent on 
the endpoints positions as well as on t\ — to. 

In what follows we will restrict ourself to the case when 
the power A is equal to 1. This simplified setup enables 
us to compare the stringy approach with the Caldirola- 
Kanai, Bateman and Kostin approaches [33 . 



A. Stringy versus Caldirola-Kanai approach 

The dissipative system under consideration is weakly 
non-Lagrangian. This means that dfi ^ 0, but fi pos- 
sesses a local integrator / ^ 0. Let us define the auxiliary 
Lagrangian function (one of plenty) as 



L CK (q,t) = e Kt \q 2 -U(q) 



(16) 



This Lagrangian is usually called CK-Lagrangian after its 
inventors P. Caldirola and E. Kanai [21 ■ Henceforward a 
damped free particle is considered only, i.e. the potential 
energy U(q) is supposed to be zero. It is immediately 
clear that the transition amplitude computed from (151 
differs from the CK-amplitude 



(<?i,*iko,*o) 



1 
1 



[2)g]exp{^ fL CK {q,t)dt} 



it, 



N 



CK 



exp{| S C K(qi,ti\q ,t )} ■ 



Indeed, the stringy and CK actions which enter the cor- 
responding transition amplitudes are different, 



'id 



4 (91 



go J 



(g + 3gi)e- Ktl -(qi + 3q )i 



-Kto 



ScA' = 



« (qi - qof 



2 e- 



G 



Since both actions are quadratic functions of the end- 
points qo and qi, and since we require that the total 
probability is conserved, the preexponential factors for 
both transition amplitudes can be computed from the 
Van Vleck formula 



1 

N 



' i d 2 S 
2irh dqi dqo 



Using this formula we arrive at the normalized transition 
amplitudes 



A lstr (q ll t 1 \q Q ,t Q ) = 



A-CK{qi,ti\q ,t Q ) = 



Airih tanh — to 



e fi 



Sic 



, (17) 



ei SoK , (18) 



2mH ( e ~ Kt ° - e" Ktl ) 



which trivially satisfy ([2|. Moreover, one immediately 
verifies that in the frictionless limit n — >• both transition 
amplitudes coincide with the free particle propagator. A 
short inspection shows that the stringy propagator de- 
pends on to and t\ only through t\ — to, but this is not 
the case for the CK propagator. The same observation 
holds for Si c i and Sck- This is a typical feature for all 
auxiliary actions defined in terms of integrators of f2, see 
|21| . We pay attention to this point because the classical 
equation of motion q — —nq which we have started with 
is invariant with respect to time translations, and one 
would naturally expect the same invariance on the quan- 
tum level. Since only Ai str (qi,t\\qo,to) has this feature, 
we obtain an efficient argument for the stringy functional 
integral ^ when compared to the CK one. The evolu- 
tion of Gaussian wave packets 

*o(«) oc e- q2 + ^ ~> cx fdq'A(q, t x | </, toxoid) ■ 
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FIG. 3: Expectation values of position (solid color lines) and 
momentum (dashed color lines) as functions of time for the 
Gaussian wave packet (the best fit for the moving classical 
particle), drawn for the stringy and CK propagators (171 
and (181, respectively. Initial wave packet characteristics 
are (q) (to) = 0m and (p)(to) = 5 ms -1 , and the friction 
constant k is set to 0.6 s . One realizes that at the time 
ii — to = In 3 ft - the stringy mean momentum becomes zero. 
Classically, for t\ — to > k _1 the physical relevancy of the 
solution breaks down. This explains the time bound on the 



applicability of stringy propagator (17 1 



rather than damped; thus, they evolve according to the 
time reversed dynamical equations. In our case we have: 

q = —Kq (damped) AND Q = +nQ (amplified) . 

These equations of motion can be derived from the least 
action principle with the quadratic and time independent 
Bateman (-Morse- Feshbach) Lagrangian: 



LB(q,q,Q, Q) = qQ + 7;(qQ - Qq) ■ 



(19) 



whose dynamics is governed by (17 1 and (18 1, is visual- 
ized in Figure |3j For more detail, see |22j . 

From the figure we can see that the stringy averaged 
momentum (p) becomes negative (i.e. meaningless) for 
ti — to > In 3 k~ . Let us try to explain this peculiarity. 
It is well known that the relevancy of the classical solu- 
tion q c i(t) of q — —nq breaks down when the time t\ — t Q 
exceeds the relaxation time k^ 1 . Since q c i is used in the 
derivation of the formula (17 1 for the stringy propagator, 



its applicability is automatically restricted as well. After 
this restriction is taken into account, the stringy evolu- 
tion of (p) becomes more acceptable from the physical 
point of view then the evolution in the CK approach. 



The canonical as well as path integral quantization of 
the theory encounters various difficulties because non- 
normalizable states of the amplified system must be em- 
ployed |23 ; . This problem will be discussed later, when 
the treatment of the auxiliary Q-degrees of freedom in 
this approach will be described. The path integral eval- 
uation of the transition probability amplitude based on 
is straightforward [U]. The result is: 

A B {qi,Qi,ti\q ,Qo,h) 



4irH sinh ^ 



x exp 
where 



{^[Qi(qi - qoP-) + Qo(qo - qiP+)] } , (20) 



B. Stringy versus Bateman approach 

The key element of the Bateman(-Morse-Feshbach) ap- 
proach \3\ are new subsidiary degrees of freedom intro- 
duced in addition to the initial ones, which arc amplified 



2 tanh^ 



Pi 



cosh ^ 



T = t 1 -t 



Our aim, however, is to find the transition amplitude for 
the damped (sub)system only. In order to obtain it we 
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must project out the nonphysical Q-degrees of freedom. 
When using the standard formula 



, trial 



(gi,*l|t?0)*o) = J dQA B (qi,Q,ti\q a , Q,to) 



we do not reproduce the free particle propagator for n = 
as desired. To overcome this trouble the following "repair- 
ing prescription" is introduced (for more detail see [25 ): 
the amplitude for the damped (sub) system to pass from 
\qo) to |<7i) within the time T is equal to the amplitude 
for the whole system to pass from the nonphysical state 
ko)|^-) to the nonphysical state |gi)|W + ), 

Af (quhlQoM :=Amp(|9,)|*_) \qi)\9+) 

dQdQ' ¥+(Q') A B (q u Q\ h\q , Q, t ) *_(Q) , (21) 



where the non-normalizible Q-system states |^±) are 
chosen as 



*±(Q) = 



'icoshfT 



=F2ia 
nft 



exp{±-aQ 2 } 



Substituting (20 1 into (21 1 we obtain the following effec- 



tive propagator for the damped (sub)system: 



A|f(g<i,ti|<Zo,*o) 



47rifitanhfT 



x exp< 



{^f [(9i-?o/3-) 2 + (go-9i/3 + ) 2 ]}. (22) 



The form of the auxiliary states |\&+) and |$_) guar- 
antees that the effective Bateman propagator for the 
damped (sub)system satisfies the normalization condi- 
tion 



The Bateman effective propagator (22 1 depends only 
on the time difference T, which favors it in compari- 



son with the CK-propagator (18 1. On the other hand 



the artificial "repairing procedure" involving nonphysical 
states discredits the Bateman approach with respect to 
the stringy one. The Gaussian wave packet character- 
istics obtained from (17) and (22 1 are drawn in Figure 
[4] As seen from the figure, the stringy quantization pos- 
sesses again better qualitative features than the approach 
to which we have compared it. 




Stringy 



FIG. 4: Expectation values of the position and momentum 
as functions of time for the Gaussian wave packet, drawn for 
the stringy and effective Bateman propagators (17 1 and (22 1, 



respectively. The initial wave packet characteristics and the 
meaning of solid and dashed color lines are the same as in 
Figure [3] 



equivalent description of the system by the Schrodinger 
equation [4]: 



ihd t ^{q,t) = 



where K[i>] is the ^-dependent Kostin potential defined 
as 



K[V] 



K 

i — 
2 



In 



*(g,*) 
*(<?,*) 



This is known as Kostin-Schrodinger(-Langevin) equa- 
tion. In the Kostin's paper a complete set of solu- 
tions of this equation in the special case of free particle 
(U(q) = 0) is presented. The set consists of the states 
ty po (q,t) labeled by the continuous quantum number po, 
the initial momentum of the particle. The state \& po (q, t) 
starts as a momentum eigenstate with the momentum 
Po, and remains the momentum eigenstate also later, but 
with the decreasing momentum, 

py po (q,t) =Poe- K(t - to) * Po (g,i). 



C. Stringy versus Kostin approach 

A possible incorporation of the classical dynamics into 
the quantum one can be obtained by using Heisenberg 
equations for the position and momentum operators |34j : 

q = p, p = F(q,p), 

where q and p obey [q,p] = iH for all t. In the special 
case when F(q,p) = —dU(q)/dq — Kp, Kostin found an 



From the ensemble of non-localized Kostin's states one 
would like to form normalized wave packet solutions, 
localized in configuration as well as momentum space. 
This is, however, impossible since the Kostin-Schrodingcr 
equation is nonlinear; thus, the Kostin approach belongs 
to nonlinear quantum mechanics and lacks superposition 
principle. In Figure[5]we plot the expectation value of the 
momentum for the Kostin state \& Po (g, t{) and the stringy 
evolved state = / «V Ai str (q,ti\q' > t )^f Po (q' \t ) 

for various values of K [35]. As k tends to zero, the stringy 
and Kostin results match each other. 
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0.1s" 1 
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0.3s" 1 
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FIG. 5: Expectation value of the momentum as a function 
of time for an eigenstate of initial momentum ty vn [a , to) oc 
exp{|gpo}, drawn for the stringy propagator (171 (blue 
dashed lines) and the propagator given by the Kostin- 
Schrodinger equation (231 (red dashed lines), respectively. 



The values of the friction constant k are indicated next to 
the curves. The initial momentum is po = 5 ms -1 . 



VI. QUANTIZATION OF THE DOUGLAS 
SYSTEM 



After we have dealt with the weak non-Lagrangeanity, 
let us consider the simplest strongly non-Lagrangian sys- 
tem. It was proposed by Jesse Douglas (one of the two 
winners of the first Fields Medals awarded in 1936) when 
studying the inverse problem of variational calculus [5]. 
The system is governed by the following dimensionless 
dynamical equations: 



-y 



(24) 



Px j Px Py j 

Py = y> fiy = -y- 

A quick calculation shows that the associated two-form 



£1 = d(p x dx + Pydy 



P, 



dt 



P 



2-*-T*) 



Pydx A dt 



does not possess a non-trivial local integrator /. Con- 
sequently, equations (24 1 cannot be obtained as Euler- 
Lagrange (Hamilton) equations. 

From the standard point of view, this curious situation 
is stalemate: no Lagrangian => no Quantum Mechan- 
ics. However, the surface integral method offers a way 
to overcome this deadlock. Using the stringy functional 
integration in the extended phase space we get 

[BE] exp-j^ f ft} oc 
s 

ti 

{~\ j (h ±2 d + Wcl ~ Wd ~ XclPycl)dt}x 
to 

ti 

[Dx Dy] exp J(\x 2 + \y 2 - \y 2 - x Pycl ) dt}. (25) 

to 



oc exp 



Here, obviously, (x c i(t),p xc i(t)) and {y c i(t),p yc i(t)) stand 
for the classical solutions of p4| matching the initial and 
final endpoints for which the quantum transition proba- 
bility amplitude is sought for, x c i{t ) = x , y c i(t Q ) = y 
and x c i(ti) = x\, y c i{ti) = yi- In the classical limit 
h — > we want the stringy propagator ([£]) to maintain 
the properties of the original dynamical system. This 
guides us to define: 



(xi,yi,tx \x , yo, to) = l(\x 2 c i+\y 2 i-\yli-x c ip yc i)dt . 



The exponential exp{jrS]j c i} of the classical action de 



fined in such a way cancels the prefactor in ( 25 



we arrive at the following transition amplituc 
Douglas system: 



Finally 
e for the 



Az)(xi,yi,t 1 |a;o,2/o,io) 



1 



y 2 - \y 2 - xp ycl )dt}.(2&) 



The path integral in the configuration space we have con- 
structed is quadratic and hence can be computed explic- 
itly. The normalized transition amplitude is 



A D (a;i,j/i,ti|xo,yo,*o) = 



2nih 



d 2 Si>ci d 2 Sud 



where 

(xi - x ) 
SDci = —. — + tan 



dx i 8xq dxi dyo 

d 2 &Dci d 2 Sjj c i 
dyx dx dyi dy 

2 I h-t \ (yi + y ) 2 



2(h-t ) 

^ (^) 

(ti - *o) 
3 

4 sin (<i — to) 
(ti - h) 



2 ) 2(h-t ) 
(xi - x )(y 1 + y ) ~ (xxyx - x y ) 

[{yi + Vo) cos (ti - t ) - 2 yi y ] 

[(y 2 + yl)-2yiyocos (h-t )] . 



4shr (ti - to) 



VII. CONCLUSION, DISCUSSION AND 
OUTLOOK 



In the paper we have developed a new quantization 
method that generalizes the conventional path integral 
approach. Throughout the paper we considered only the 
nonrelativistic quantum mechanics of spinless systems. 
However, the generalization to the field theory is rather 
straightforward. We have just to pass from the space of 
particle positions to the space of field configurations. 

The mathematical language of our exposition respects 
the Vladimir Arnold "principle of minimal generality." 
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Formulas are mostly written down in one (local) chart. 
From this, however, one can ascend to a global, coor- 
dinate free description employing bundles and jet pro- 
longations. In order that we did not cloud up the main 
idea of the paper, we did not follow such fluffy approach. 
However, in the future it could be useful to analyze ob- 
structions to the implementation of our method which 
arise from the global properties of the underlying space- 
time geometry. 

Special attention was paid to the stringy quantiza- 
tion of dissipative systems and of the Douglas system. 
The latter cannot be quantized by any known quantiza- 
tion technique while the former can. In fact, dissipative 
(friction force) systems were studied extensively in the 
past and several approaches to their quantization were 
proposed. The stringy quantization was compared with 
three of them, Caldirola-Kanai, Bateman and Kostin, 
and it was shown that it either gives better results (in 
the first two cases) or is simpler to apply (in the third 
case). 



B. Surface functional integral versus 
Feynman- Vernon 

The goal of our paper was to write down a functional 
integral formula for the quantum transition amplitude 
in terms of the underlying classical equation of motion. 
However, we did not eliminate the artificial notions of 
pure states and classical action for non-Lagrangian sys- 
tems. In the surface functional integral (|9J, these notions 
were implicitly present. 

A possible alternative to our approach consists in ex- 
pressing the transition probability (not the probability 
amplitude!) using the functional integral. The prob- 
ability that the system evolves from the mixed state 
Poilo-i q'oito) at the time to to the mixed state p ± (qi, q[,ti) 
at some later time t\ is 

Pp„^ Pl <* dq d(/ dqidq[p (qo,q' ,t )p 1 (qi,q'i,ti)>< 



D7][Df]exp{^(7 - f){ P dq-Hdt)} 



A. Stringy versus Caldeira-Leggett model 

As mentioned in the introduction, the quantization 
proposed here shares a common property with the well 
established approaches of Caldirola-Kanai, Bateman, 
Kostin and others: they all fail to cover decoherence phe- 
nomena. On the other hand, as we have seen, the stringy 
quantization produces a quantum propagator that does 
not respect the memoryless Chapman-Kolmogorov equa- 
tion. We believe that this essential feature of the theory 
is measuring/reflecting the non-Lagrangeanity of the sys- 
tem on the quantum level. 

From the author's point of view one can rise a concep- 
tual objection also against the particle-plus-environmcnt 
model by Caldeira and Leggett [5]. Their approach is 
microscopical except the way in which the spectral den- 
sity pd(u) for the reservoir degrees of freedom is intro- 
duced. This density is not obtained from any kind of 
microscopical theory and the form it assumes is moti- 
vated solely by the necessity to obtain the dissipative 
term —nq in the effective theory. On the top of it, the 
interaction Hamiltonian is such that it produces, after 
integrating out the reservoir degrees of freedom, an ad- 
ditional Langevin stochastic force F(t). Thus, the effec- 
tive classical motion of the particle is governed by the 
Brownian equation of motion p = —up + F(t), which is 
conceptually different from the physical situation consid- 
ered in section [V] It is also worth to point out that the 
Caldeira-Leggett model cannot describe satisfactorily the 
friction force proportional to — p A with 1. There is 
no doubt that for A = 1 the model gives better results 
than all effective theories mentioned before including the 
stringy one, but it is fair to say that it describes a differ- 
ent physical phenomenon, namely the quantum Brownian 
motion of a particle which is in thermal equilibrium with 
a reservoir. 



Here 7 and 7' are curves in the extended phase space 
whose g-projections connect qo with q\ and q with q[ 
respectively. Using the Stokes theorem one is able to 
convert the difference of the line integrals of the one- 
form pdq — Hdt into the surface integral of the two-form 
Q = d(pdq — Hdt). Explicitly, 



- J )pdq-Hdt=- J Q-U - J jpdq, (27) 

77' E Ai Ao 

where E represents again a map from the parametric 
space (t, s) € [to, ti] x [0, 1] to the extended phase space, 

£ : (t, s) h-> E(f, s) = [qs(t, s),ps(f, s),ts(t, s) =t) , 

such that <9E = •y' (t) — Xi(s) — j(t) + Xo(s) . The sideways 
A-boundary curves [36] of £ 

Ao(s) = £(io, s ) an d \%(s) = Yl(tx, s) 

live in the instant phase spaces of M. 2n+1 [q a ,p a , t], i.e. 
in the submanifolds R 2n [q a ,p a ,t Q } and R 2n [q a ,p a ,ti\. 
Moreover, as the intrinsic parameter s varies from to 
1, the gs-components of Ao(s) and Ai(s) vary from qo to 
q'o and from qi to q[ respectively. 

Let us denote by U the space of all £-maps for which 

ga(*o,0) = q , qs(h,0) = qu A (s) C R 2n [q,p, to], 
<7 S (io,l) =<?(,, <te(*i,l) = «i, Ai(a)cR 2n [g,p,ti]. 



Since the right hand side of ( 27 1 depends only on 



E and its two boundaries, the double path integral 
J [D7] [D^f 1 ] ■ ■ ■ entering P (9o _ >Pi can be rewritten as a 
surface functional integral J[DE] • • • , 



p ->-Pi 



x oc ldqodq'odqidq' l p {qo,qo,to)p 1 {qi,q' x ,tx) x 



DE] exp{~[/"r!+ ( f-hpdq\y (28) 
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The formula we have just arrived at was derived under 
the assumption that f2 = d(pdq — Hdt). However, it is 
clear that the expression in the exponent is free of any 
reference to the Hamiltonian and requires just the clas- 
sical equations of motion. Therefore it seems reasonable 
to postulate the probability formula ( 28 1 also for non- 
Lagrangian systems. 

The approach we shortly presented here does not use 
the notion of pure states of a non-Lagrangian system. It 
resembles the Fcynman- Vernon approach [10], in which 
one introduces the influence functional in the presence 



of dissipative forces. Our probability formula (28) uses 
the surface functional integral in the extended phase 
space, while in the Feynman- Vernon approach one just 
has to calculate a double path integral in the configura- 
tion space. There is a chance that after a certain dis- 
cretization we will be able to convert the surface func- 
tional integral into the double path integral in the con- 
figuration space (see appendix|B|, and as a result, we will 
find the explicit form of the influence functional. Work 
on this topic is in progress. 
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(9»,to) 



FIG. 6: Vector field W acting infinitesimally in the extended 
phase space. As a result, the initial surface E ea;t G U^^) 



G U, 



and its bound- 



A. 



M. 



D. 



moves to some nearby surface T, ext c - < - - j , 
ary dH ext = 71 — Ai — 70 + A0 is transported to a new boundary 
<9Ef^ = ! y' 1 — A'j — 7q + A|j. Since the variational vector field is 
assumed to preserve II, W stays tangential to the momentum 
submanifolds with fixed (qo,to) and (qi,ti). 



IF-fields it holds 

o = iim §(s ^ ) ~ s(Sezt) = /r u o. 

<5->0 5 



Appendix A: Stringy variational principle 

In section |III| w e shortly outlined how the transi- 
tion amplitude (|9j can be related to the stringy varia- 
tional principle, with the action S defined in (111. In 



what follows we will perform a complete variation of S, 
with the domain extended from U 7ci to a wider class 



U 



Ik U, u 



(71,70)- 



The stringy set U contains all 



surfaces in tie extended phase space trapped between 
the submanifolds with fixed (go,^o) an d (<7i,ii) (includ- 
ing separate histories, regarded as shrunk surfaces). In 
the course of variation the classical dynamics (|5| we have 
started from will be recovered. 

Suppose we have an extremal surface Y> ext G U( 7l C 
XL and a variational vector field W defined in its neigh- 
borhood such that the flow of W preserves U. By defini- 



tion, W moves Yi„ xt G 



(71,70) 



to some other worldsheet 



SfS' G U,-' -'x, where 8 is an infinitesimal increment of 

ext (7i,7 ) 

the parameter of the flow generated by W. The situation 
is depicted schematically in Figure [6j 

The extremality of Y, ext means that for all variational 



where Cw stands for the Lie derivative. Let us compute 
the surface integral on the left hand side: 



J c w n = J d(Wj n) + Wj dn 



= JWjn+JWjdfl, (Al) 

where the the symbol j denotes the inner product (con- 
traction) of a vector with a differential form. The first 
term on the right hand side depends only the values of f2 
on the boundary of T, ext and can be recast into the form 



The last two terms here give individually zero contribu- 
tions, as can be seen from the following quick consider- 
ation: the vector field W is assumed to preserve U and 
therefore its restrictions to the A-boundaries are 



w \\= W a{P><lo,toj .. 

IA \ > Q 



d 



w\ x =w a ( P>qi ,h) 



d 
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The one-form we are integrating is W j — W a (dq a — 
(p a /m)dt), but since both q and t stay unchanged on A. 
the integral is zero. As a result, to annihilate the bound- 
ary term for all variational fields W we are forced to chose 

the curves 7q = (£fo(*) ? Po (*)>*) and 7i = (&(*), Pi (*),*) in 
the extended phase space in such a way that their instant 
tangent vectors annihilate f2 at any moment t G [to,ii]- 
Thus, it must hold 



(A2a) 



or, equivalently, 



0= (F a -p a )dq a 



q a --)dp a + 
m J 

q a F a - Pa ^)dt, (A2b) 



for both 70 and 71. Here one recognizes the dynamical 
equations (|5l as desired. 

From the boundary term analysis and the assumption 
about the uniqueness of the history 7^ which we adopted 
from the very beginning we can conclude that E e;r4 G 



C lt 7ci . Hence E e;Et has the topology of a closed 



(7c!,7cl) 

string attached to the classical trajectory 7^. To specify 
its shape the second term in ( Al I must be employed. The 



resulting variational equation, written in a coordinate- 
free notation, is 



= dil (ftEeat, <9 s E e 



(A3a) 



This is equivalent to the following system of partial dif- 
ferential equations for the unknown functions q a {t, s) and 
p a (t, s) (indices a and k run from 1, . . . , n): 







= 



dq k dF k 
ds dp a ' 

V fdFa _ dFk 
ds \ dq k dq a 



dpk dF a 
ds dp k 



(A3b) 



One solution of these equations, satisfying all boundary 
conditions, is trivial. It is the shrunk surface E ex t(t, s) — 
"f c i{t)- After recalling the assumption of uniqueness of 
7ci once more, we can see that this is the only solution of 
( A3 ) . If there existed a closed unshrunk extremal surface 
Eeici, the initial set of equations of motion (|5| would have 
at least one one-parametrical family of classical solutions 
between the given pair of endpoints. 

The variational principle we described above operates 
on a wider stringy class IX than it is in fact necessary. The 
surface integral formula ^ requires just the restricted 
subset U 7o! C It. The transition to U 7cl is advisable for 
two reasons. First, in the Lagrangian case with 57 = 
d(pdq — Hdt) we obtain an equivalence (modulo additive 
constant) between the least action principle using § and 
the standard Hamilton least action principle. Explicitly, 

E G l% 7ci) c U 7c! i y S(E) = (/"/) (P d i - Hdt ) ■ 

7 id 



Second, the wider class It contains plenty of degenerate 
(shrunk) surfaces, the histories 7. These are obviously 
stationary surfaces of § when being varied within It, since 



Hi 5w ) - S(7) 



n 



Q=0-0=0 



for any variational vector field W. However, only one 
of these histories, namely 7,^, satisfies both equations 
(A2| and (A3) at the same time. This fictitious problem 



is avoided when one works from the beginning with the 
stringy subclass U 7ci C II. 



Appendix B: Surface functional integral — 
computational details 

Let us explain here in some detail how the surface 
functional integral is computed. To be as tangible as 
possible consider a one-dimensional system only, so that 
the extended phase space will be the three dimensional 
space M. 3 [q,p, t] (more dimensions represent only a tech- 
nical problem). The particle is supposed to move under 
the combined action of the potential U(q) and the fric- 
tion force with a general A-power law. The dynamical 
equations are 



q=p, p 



U{q) 
aq 



up 



and in addition to them, we require that the particle 
satisfies the boundary conditions q(to) = qo and q[t\) = 
qx . According to the definition Q , the two- form O is 

Q = d(pdq - l,P 2 dt - U(q)dt) - np A dq A dt . (Bl) 
Our aim is to compute the surface functional integral 



9l over the stringy set U 7cl = (J - It 



J 7 u (7,7 ci )' 



The direct 



application of (10 1 and ( B 1 1 together with the Stokes 
theorem yields 



(B2) 



[DS]exp{^ j n} 

s 

[D 7 ]exp{^( 
DS] exp{— j 



(pdq- ~p 2 dt-U{q)dt)} 
p A dq A dt j . 



The nontrivial part of this is the functional integral over 
the stringy subset lt( 7j7ei ). As mentioned earlier, S € 
U( 7j7cl ) is a map 



E : [io.ii] x [0,1] -> 
such that for Vi € [to, ii 
E(t,0) =7ci(t) 



R 3 [q,P,t], (i,s)^E(i,s) 
and Vs € [0, 1] there holds 



E(i,l)=7(i) 



and 



E(i ,s) G 
E(i x ,s) G 



i 3 [qo,p,t ] , u „. 
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FIG. 7: Rectangular nodal web in the parametric space 
[to, ii] x [0,1] of a surface E 6 Uf~~ ,). Each elementary 
tile encloses the area A ■ e = (ti — to)/(KL) (at the end, the 
numbers K and L will be sent to infinity) . The points marked 
by crosses are constrained by the conditions ( B3 1 . 



To proceed further with the functional integral in ques- 
tion, introduce a set of regularly distributed nodal points 
in the parameter space, 

{(t +rA,as) €[t ,h] x [0,1]}. 

The points of the set are labeled by two discrete indices, 
the time index r = 0, . . . , K and the space index a — 
0, . . . L, see Figure [7] In this way we obtain for any E- 
map from VL^ tlc[ ) a discretized (K + 1) x (L + l)-tuplc 



(K,L) 
(r,ff)=(0,0) 



The boundary values are required to be consistent with 
( B3 1 . Explicitly, for all indices r and a it must hold 



Q(t,o) = Qci{to + tA), <7 (7%i) = q(t + tA), g (0 , CT ) = q , 
P(t,o) = Pd(to + tA), P( r ,V) = p(t + tA), q( K ,a) = Qi- 

The discretization of E enables us to approximate the 
integral of p A dq A dt as follows [37_: 

pHq A * := £ g As {pf T , ff) } = 

T=0 <T = 

K-l L-l 



= EE 4 {?(r,a)(P(U-l) -P{r,«))) 
K-l 

+ Yl A {l{r,L)Ptr,L-l) ~ <Z(r,0) P(t, ) } ■ 



(B4) 



Formally, the discretized functional integral over all 
stringy configurations from It(^ 7c( ) is a multiple integral 
over all unconstrained variables £(t,o-),P(t,<t) which are 
needed to specify E. The only problematic part, as usual, 
is the choice of an appropriate integration measure. For 



the reason that will become clear in a moment, we choose 
the measure as 



-foo 



K-l L-l 



The first step, when dealing with the discretized func- 
tional integral 



Int= i [DE]cxp j-^K-f p A dqAdt} 



consists in the integration over the internal stringy posi- 
tions dq( T rT y After this integration, the integrand trans- 
forms into a chain of delta functions, 



K-l L-l 



int = r n n ( A pV)) d p^) s «^)-pU) = 

-oo T=1 a=1 

. K-l 

X eXP {~^ K 5I A hr,L)P(r.L-l) ~ q(Tfi)P( T ,0)]} 



The next step is the integration over the internal stringy 
momenta dp^^y After performing this trivial integra- 
tion we arrive at the expression 



K-l 



Int = cxp A[<1(T,L)P( T , ) ~ ?(r,0) Pfrfi)} } ■ 

r=l 

In the exponent there appears the discretized version of 
the integral 

J dt{q(t)pii(t)-qci(t) P ^(t)}, 

to 

where q{t) stands for the g-projection of ^{t). Conse- 
quently, after returning back to the continuum limit we 
obtain the following important result: 

[DE] exp J P A dq A dt} = 

s 

ti 

exp{- l - K [ dt[q(t)pi(t)-q cl {t)pi{t)]y (B6) 



After substituting (B6) into the initial formula (B2) we 



recover a path integral in the extended phase space. The 
integral is quadratic in momenta, with the discretized 
(standard) Liouville measure 



K-l 



dpK tt dpr dq T 

T=l 
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The integration over momenta can be carried out explic- 
itly and we finally obtain a path integral in the configu- 
ration space only (tildes are removed from the position 
variables) , 



classical action ( 13 1 



[2)E]exp{^ j n} = 

s 

ti 

cxp[-^ J [^q%-U(qci) - Kq c ip£)dtj X 

to 

tl 

[D<Z]exp{| J {^e-U{q)-Kqpi)dt). (B7) 



This formula served us as a guide when introducing the 



Let us make one final comment. In our functional 



measure (B5) no integration over the momenta P(o,cr) 
and P(k,<j) was prescribed. These momenta are how- 
ever needed when the boundary of the discretized surface 
£ G U 7ci is specified. They define the auxiliary curves 



Ao(s) and Ai(s) introduced in section III The curves were 
chosen completely arbitrarily, but as seen from ( B4 1 , the 
result is not affected by them at all. Hence, since every- 
thing substantial is independent of P(o,cr) an d P(K,a)- it is 
justifiable to discard these quantities from the functional 
measure [f£]. If we would not do that for some reason, 
they will integrate into an artificial infinite factor, which 
will be removed anyway after we apply the normalization 
conditions (|2j. 
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~ K QPci—KL (t) and 



-nqPd-KM 



KL 

(t) for t € [t K ,t L ] 

LM (^) for 1 £ [tL,t M ] 
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